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The evolution and viscous damping of cosmic magnetic fields in the early uni- 



D ■ verse, is analysed. Using the fact that the fluid, electromagnetic, and shear 

W . viscous energy- momentum tensors are all conformally invariant, the evolution 

^O , is transformed from the expanding universe setting into that in flat spacetime. 

Particular attention is paid to the evolution of nonlinear Alfven modes. For 
^ I a small enough magnetic field, which satisfies our observational constraints, 

these wave modes either oscillate negligibly or, when they do oscillate, be- 



f^ . come overdamped. Hence they do not suffer Silk damping on galactic and 

^^ , subgalactic scales. The smallest scale which survives damping depends on 

f~^ \ the field strength and is of order a dimensionless Alfven velocity times the 

^^ \ usual baryon-photon Silk damping scale. After recombination, nonlinear ef- 

fects can convert the Alfven mode into compressional, gravitationally unstable 
waves and seed cosmic structures if the cosmic magnetic field is sufficiently 
strong. 



I. INTRODUCTION 

The origin of ordered, large-scale cosmic magnetic fields remains a challenging problem. 
It is widely assumed that magnetic fields in astronomical objects, like galaxies, grew by 
turbulent dynamo amplification of small seed magnetic fields [|T|. There have been many 
suggestions for producing the small seed magnetic fields required to prime the amplification 
process [^ . Many of these proposals appeal to processes which might operate in the very 
early universe . Moreover, the efficiency with which known turbulent dynamo mechanisms 
can produce the observed fields is still being debated 0]. An alternate possibility is that 
the galactic field could be the remnant of a larger cosmological field of primordial origin 
^, although, as yet, there is no compelling mechanism for producing the required field. 
It could form part of the initial conditions, arise at a phase transition, or be produced in 
some way at the end of a period of infiation [Q]. If the primordial field is homogeneous then 
the isotropy of the microwave background places a limit of 3.4 x 10~^(r2o/i5o)2 Gauss on its 
present strength . This limits the strength of any primordial field on scales greater than 
the horizon at last scattering. 

Dissipative processes play a key role in fashioning the spectrum of irregularities that 
survive the radiation era. It is well known that density fiuctuations in the baryons suf- 



fer "Silk" damping in the early universe, due to radiative viscosity 0. This leads one to 
wonder whether fluctuating magnetic fields produced in the early universe, which couple 
to the baryons, will also be damped by similar processes. This problem was considered by 
Jedamzik, Katalinic and Olinto P], who examined the damping of linearised magnetohydro- 
dynamical (MHD) wave modes in an expanding universe. However, there remain subtleties 
about the formulation of the problem, and the damping processes, which prompt us to re- 
analyse the problem in a completely different fashion that also allows non-linear effects to 
be considered. 

We are able to do this by exploiting the conformal invariance of the relativistic fluid, 
electromagnetic, and shear viscous energy-momentum tensors to transform the problem from 
the expanding universe setting to a simpler problem in flat spacetime. This transformation 
is explained in detail in the next section and in Appendix A, and it allows us to study the 
evolution and damping of various MHD waves. After taking the non-relativistic limit of the 
resulting MHD equations, in Section HI, we focus in Section IV on non-linear Alfven wave 
solutions and study their damping. The compressible MHD modes are also briefly treated, 
in their linearised limit, in Section V. As the universe expands, the relevant mean-free-path 
of the particle responsible for the damping eventually exceeds the wavelength of a given 
mode. The evolution of various MHD motions in this free-streaming regime, is considered 
in Section VI. Eventually, the universe recombines to form atoms and the tangles in the 
magnetic fleld that survive damping, can create gravitationally unstable perturbations. This 
post-recombination regime is considered in Section VII. In Section VIII, we connect these 
ideas together to discuss how a given spectrum of magnetic inhomogeneities evolves, and 
section IX summarizes our conclusions. 

II. CONFORMAL INVARIANCE OF RELATIVISTIC, VISCOUS MHD 

We shall make use of a conformal invariance property of the fleld equations of general 
relativity. Consider two metrics which are related by conformal transformations of the form 

9% = ^'9,u. (2.1) 

where Q = Q{x"). Suppose the cosmic fluid is described by an energy-momentum tensor, 
say T^'^. If this is traceless, Tj^ = 0, then the equations of motion of the fluid are invariant 
under conformal transformations. That is, we have 

Tf"".^ = ^ T*^;',^ = with T*f"' = Q-^T^"". (2.2) 

Here ; * denotes a covariant derivative with respect to the metric g*. This result is straight- 
forward to demonstrate by direct calculation. 

We apply this result to the early universe setting. Firstly, note that the Friedmann- 
Robert son- Walker (FRW) metric is conformally flat. We will only need to consider the 
zero-curvature metric with line element 

ds"^ = g^ydx^dx'^ = a^ir) —dr"^ + dx^ + dif + dz^ , (2.3) 

where r is conformal time and a{T) is the scale factor. We also deflne the comoving proper 
time t using dt = adr. We can relate this metric to flat spacetime, r/^j,, by a conformal 



transformation, ri^i, = (7*^ = a~^{T)g^i,, with fl = a~^{r). So if the energy-momentum 
tensor of the matter in the early universe is trace-free and transforms under conformal 
transformations as in Eq. ( p.2| ), then we can transform the fluid equation of motion to 



a simpler equation in flat spacetime. This approximation method allows us to transform 
non-linear solutions to the MHD equations, for a fluid with trace-free energy-momentum 
tensor, from flat spacetime into an expanding FRW universe. We are, of course, neglecting 
the back-reaction of those motions upon the form of the metric. That is, the gravitational 
potential (metric) perturbations due to the inhomogeneous motions must be small for this 
approximation to hold. For example, if the inhomogeneous motions of a radiation fluid give 
rise to energy density inhomogeneities of amplitude Sp/p over a length scale A, then the 
metric perturbations induced are of order {6p/p){\/ct)^. These are small because, for the 
parameters we consider, Sp/p << 1 and the perturbations are on sub-horizon scales. 

The energy-momentum tensor in the early universe will be modelled by a linear combi- 
nation of ideal fluid, non-ideal fluid, and electromagnetic fields. We have 

rpt^U ^ ^M. ^ J^U.^ ^ yM^^ ^2.4) 

where we have separated the energy-momentum tensor into an ideal fluid part Tj''^ , a non- 
ideal fluid part T^j and an electromagnetic part T^^. We take the ideal part to be a perfect 
fluid with equation of state p = p/3, where p is the pressure and p is the energy density 
(we have set the speed of light c = 1.) For most of the period before decoupling, the early 
universe is radiation dominated and one can use the above equation of state. Later, we will 
comment on the effects of matter-domination {p = 0) before decoupling. Thus, we have 

T^ = ip + p)^^ + pg"". (2.5) 

where If^ is the normalised four-velocity of the plasma, satisfying the condition U^U^ = —1. 
Note that for p = p/3, we have Tj" = 0. Under a conformal transformation, the components 
of Tj''^ transform to a new set of (starred) variables as follows 

p* = Q-^p^ p* = fi-4p, t/*M = Q-ijjf^ (2.6) 

and the ideal fluid energy- momentum tensor transforms as Tj'^'^ = Vt~^Tj'^ . 

The non-ideal fluid part of the energy-momentum tensor can be written as 

TZ = -vH^^'H-'^W^p, (2.7) 

where 

and 1] is the effective shear viscosity coefficient, given by (see Weinberg 0); 

where T is the temperature of the radiation-dominated universe, g the statistical weight and 
Id the mean-free-path of the diffusing particle. In the radiation-dominated epoch, the bulk 



viscosity is zero and we have neglected the thermal conductivity term since it does not affect 
the Alfven wave mode that we will focus upon. 

One can easily check that the trace T^j = 0. Also, under conformal transformation, 
using U*'' = n-^W and U*^ 
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and H*^""!!*^^ = Q 4^/ia^;//3. g^ ^^le non-ideal stress conformally transforms as 
TX = -v*^'^H^"'H''^W^p = Q-^TZ, if V* = ^'^V- 
The electromagnetic part of the energy-momentum tensor is given by 
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Pl^pyi _ -gf'^F^sF^^ 



(2.10) 
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where F^^ = A^.^ — A^.y = A^^^ — A^^, is the Maxwell tensor, A^ the four-potential, and a 
comma denotes an ordinary derivative. This is also traceless: T'^^j = 0. The evolution of 
the electromagnetic field is governed by the Maxwell equations 



F"".,, = 47rJ^, 
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(2.13) 



Here, J^ is the four-current density. Under a conformal transformation, these equations are 
also invariant provided we transform the electromagnetic field tensor and current density as 
follows: 



Under these transformations we verify that 
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Finally, consider the relativistic generalisation of Ohm's law. 



(2.14) 



(2.15) 



(2.16) 



where a is the conductivity of the fiuid, measured in the fiuid rest frame. Transforming 
the four-current and velocity to the starred frame, we can see that Ohm's law also remains 
invariant under conformal transformations if we define a starred conductivity, a* = Q~^a. 

/^From these preliminaries, we see that the total energy-momentum tensor consisting of 
an ideal fiuid, non-ideal fiuid, and electromagnetic field is trace-free, and also transforms as 
in Eq. ( |2.2| ) under a conformal transformation. Suppose we choose Q = a'^lr), with the 
scale factor normalised such that a(ro) = 1 at present conformal time Tq. Define a new set of 
starred variables as in Eqs. ( |2.tj| ), ( |2.11| ) and (2.14). The evolution equations of the starred 
variables then become, 
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(2.17) 



Since these equations are simply the flat spacetime equations of relativistic, viscous MHD, we 
can carry over all results obtained in the flat spacetime context to MHD in FRW spacetime 
on scales such that the metric perturbations created by the magnetic inhomogeneities and 
their resulting motions are small. The only extra complication introduced by cosmology is 
the time-variability of the viscosity and conductivity coefficients. (The transformation of 
MHD equations from the FRW universe to flat space, has also been explicitly constructed 
in ref. ( [|T3])- However, they did not include the viscous stress, nor exploit the simple 
conformal invariance properties stressed here. The idea that flat space solutions can be 
taken over to FRW universe in the radiation era, has also been used by Liang |]TT| , to study 
shock waves in the radiation era; but, without the inclusion of magnetic flelds). 

At this stage, it is useful to point out that any particle species which does not interact 
with the relativistic fluid is also easy to take into account. Such a non-interacting component 
may be any species of dark matter, or neutrinos, after they have decoupled from the rest 
of the matter. The energy-momentum tensor for this component will then be separately 
conserved. As long as this component perturbs the form of the FRW metric negligibly, the 
above treatment of conformal transformation of relativistic MHD to flat spacetime will go 
through. This also means that the MHD equations in the above form can be applied to the 
baryon-photon-magnetic fleld system, even after dark matter domination, until the baryon 
density becomes comparable to the radiation density, which happens close to decoupling. 
In the sections to follow, we will use these equations to examine the evolution of MHD 
waves, in particular the non-linear Alfven mode, after taking the non-relativistic limit of the 
equations. 

III. NON-RELATIVISTIC LIMIT 

Suppose the universe was seeded with magnetic flelds at some early epoch in its radiation- 
dominated phase. We consider the subsequent evolution. For magnetic flelds of realistic 
strengths (ie that are not going to produce highly discordant microwave background fluctu- 
ations or radically change the course of primordial nucleosynthesis) , the induced bulk fluid 
velocities are in general highly non-relativistic. Suppose we write the four-velocity in the 
transformed metric as U*^ = (7, 7v), where 7 = (1 — v^)~^/^. This form satisfles the normal- 
isation condition on the four-velocity and v = rfx/rfr = arfx/rft is exactly the peculiar bulk 
velocity of the fluid in the FRW metric. The non-relativistic limit corresponds to taking 
|v| << 1 in the fluid conservation equation in Eq. ( |2.17| ) . 

If we deflne the electric fleld E* = {E*^ , E*^ , E*^) and the magnetic fleld B* = 
(_B*\ 5*^, B*^) in the starred metric by 



then the time component of the fluid equation in Eq. ( 2.17 ), in the non-relativistic limit, is 

^ + V.[(p* + p*)v] - E*.r - r]*V.i = 0. (3.2) 

or 

Here, J* is the current-density vector whose components are the spatial components of the 
transformed four-current J*'^, and f = V(v^/2) — (2/3)vV.v. The spatial components give 
the Euler equation 



^ [{p* +P*)v] + (v.V) [(p* +p*)v] + vV. [(p* +p*)v] 

= -Vp* + r xB* + r]* 



VV + -V(V.v) 

o 



(3.3) 



We have assumed here that the net charge density (J*°) is neghgible. 
The Maxwell equations in the starred metric, are 

V X B* = AttJ* + -—- V X E* = — —- V.B* = V.E* = 47rJ*° (3.4) 

OT or 

They are supplemented by the non-relativistic limit of Ohm's law 

E* + vxB* = — . (3.5) 

a* 

The Maxwell equations in a "Lab frame", with fields E and B defined using the fundamental 
observers of FRW spacetime, are also easy to write down and are given in Appendix A. 

In what follows we shall make the usual assumption that, for non-relativistic velocities, 
the displacement current term can be neglected. In this case, we can split up the Lorentz 
force term in the canonical way as J* x B* = — VB*^/(87r) + (B*.V)(B*/47r). We also 
assume that the early universe was a perfect conductor (see for example ref |0]) and take 
the limit of a oc a* ^ oo, in Ohm's law. The magnetic field then satisfies the ideal limit of 
the induction equation 

^ = V X [v X B*] . (3.6) 

or 

In terms of the "Lab" magnetic field, defined in Appendix A, the ideal-limit induction 
equation is 



^7^ — - = — ^^ — - = -V X 
at a OT a 



V X (a^B)] . (3.7) 



Let us turn now to consider solutions of the above equations. 

Consider first the unperturbed evolution, with a zero peculiar velocity and negligible 
magnetic fields. The solution of the fluid equations is then p* = p*/3 = Ci with Ci a 
positive constant. This implies that, in the original variables, p = p/3 = Ci/a^, as expected 
for the radiation-dominated universe. (Furthermore, for zero peculiar velocity, for any "test" 
magnetic field, B* is constant in time, or B oc a~^, a result which is intuitively expected for 
the "Lab" magnetic field due to flux freezing in the expanding universe.) Now consider the 
effect of introducing tangled magnetic fields in the universe. The Lorentz force associated 
with a tangled field will cause the fluid to move and induce a non-zero peculiar velocity. The 
coupled system of equations describing the evolution of the velocity and magnetic fields is 
highly non-linear. For this reason, the authors of ref. p] examined only the case of weak 
perturbations around a quasi-uniform field. We shall to begin with, follow a complementary 
approach and look at special non-linear solutions. This will also give some feel for how a 
general tangled field configuration will evolve. 



IV. NON-LINEAR ALFVEN WAVES IN THE EARLY UNIVERSE 

A. The ideal non-viscous regime 

At sufficiently early times, (or, equivalently, for the field and velocity on sufficiently large 
scales), one may assume the matter is a perfect fluid and neglect any viscous effects. Also, 
the fluid radiation pressure in the early universe is typically much larger than the magnetic 
pressure, for the field strengths we are considering. Their ratio is given by 

= ^3x10^^5%. 4.1 

where -B_9 is the present-day magnetic field in units of 10~^ Gauss. Here, we have assumed 
that the fields are simply frozen into the uniformly expanding universe, neglecting the effects 
of the peculiar velocity. Since B*^ / {Sirp*) << 1, to an excellent approximation, one can take 
the motions induced by the magnetic field to be almost incompressible, with p* + B*"^ / {Sn) si 
p* = constant and drop the pressure gradient term in the reduced Euler equation (|3.3|) . In 
this limit we have filtered out "fast" compressible motions and we can look for solutions 
with V.v = 0. Equation ( p.2|) then gives, in the ideal and non-viscous limit, dp* /dr = 0. 
Also, the Euler equation reduces to 

<9v , „, (B*.V)B* , , 

while the induction equation is 

+ (v.V)B* = (B*.V)v. (4.3) 



A7f{p*+p*) "^r 



dr 
These equations have a wide class of exact, stable solutions of the form |T^ 

R* r)v 

0, V.v = (4.4) 

In these solutions the nonlinear terms in Eq. ( ^4.2| ) and ( [4.3|) exactly cancel. Any complicated 
tangled field pattern is possible, if accompanied by a velocity along the magnetic field at the 
local Alfven speed. 

A particular case of this solution is one where the magnetic field is split up as B* = 
Bg + b*, into a uniform and constant component, Bq, and a tangled component b*(x, r). 
Fix the co-ordinates such that Bq lies along the 2;-axis, that is Bq = BqZ, where z is the 
unit vector along z. Now choose v = b*/(47r(p* + p*)Y^'^', the nonlinear terms in the Euler 
and induction equations cancel out and we have 

dh* dh* dv dv , , 

OT oz or oz 

where we have defined the Alfven velocity by 



Here, for the numerical estimate, we have taken the dominant contribution to the energy den- 
sity p = p^, the photon energy density, as would be appropriate in the radiation-dominated 
era, after the epoch of e~^e~ annihilation. The general solution of Eq.( |4.5| ) is therefore 
a nonlinear Alfven wave travelling antiparallel to Bq, with v = b*/(47r(p* + p*)Y^'^ = 
F{x, y, z + Vat) with an arbitrary function F. One can also have another class of solutions 
with V = — b*/(47r(p* +p*))^/^, where the wave travels in the same direction as Bq. Both 
these solutions are stable and they exist as long as the two sets of waves do not overlap in 
space 
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It is not possible to generalise these nonlinear solutions to the diffusive viscous regime, in 
the above form, for arbitrary viscosity and conductivity coefficients. However, there exists 
a special case (see Landau and Lifshitz ||13|), with the velocity and tangled magnetic field 



having arbitrary strengths, but aligned perpendicular to Bq and depending only on z, where 
such a generalisation is possible. We now look at this case in more detail. 

B. Nonlinear Alfven waves in the viscous regime 



We begin by reinstating the viscous term in the Euler equation ( p.3|) . As before, assume 
that the magnetic field can be written as B* = Bq + b*, with a uniform Bq. We assume b* 
is perpendicular to Bq, but do not put any restriction on the strength of b* so that it need 
not be a sma// perturbation of Bq. We also take the peculiar velocity v to lie perpendicular 
to Bq and assume that all the variables depend only on z and r. In this case, the velocity 
perturbation automatically satisfies V.v = 0. Further, the ratio of the magnetic energy 
density to the fiuid energy density, B^/[8np) ~ lO^^-B^g << 1. So even when there is strong 
damping of motions induced by the field, and a significant fraction of the field energy density 
goes into heat, p will be perturbed negligibly. It is an excellent approximation to neglect 
the viscous term in Eq. ( p.2|) . Then, in the ideal limit, Eq. ( p.2| ) implies {dp* /dr) = 0. The 
non-linear terms in the Euler and induction equations are individually zero because there is 
no variation of b* and v along the fields. These equations then reduce to 

dr \p*+p*J L -I \An{p*+p*)l dz {p* + p*) 





dh* 


P 1 -« {01T 1 A f ^ , :.\ 

L J \47i{p*+p*)J 


dz 


dh* dv 

dr ""^ dz 





(4.8) 

Note that the LHS of Eq. ( ^.7] ) has zero divergence. The RHS will also have zero divergence 
only if the total pressure p * +B*'^/{87r) is uniform in space. As mentioned above, this is 
likely to be a good approximation for this mode since the radiation pressure in the early 
universe is typically much larger than the magnetic pressure. One can therefore drop the 
pressure gradient term in the reduced Euler equation ( [4.7|) . Writing b* = bo{T,z)n and 
V = vo{t, z)n, eliminating vq from Eqns. ( [4.7D and ( |4.8| ), gives a damped wave equation for 
&o(r,z). 



dr^ (p* + p*) dz"^ \9t I dz"^ 



where we have defined the Alfven velocity, Va, as before. This hnear equation generahses 
the nonhnear Alfven mode to the viscous regime. It can easily be solved by taking a spatial 
Fourier transform. For any mode bQ{T,z) = /(r)e**^^. we have 

f + P^f + eV^f = (4.10) 

which is the equation for a damped harmonic oscillator. 

The behaviour of solutions to the damped oscillator equation depends on the relative 
strengths of the driving and damping terms. Suppose we define 

u:, = Wa, D = ^^ocr^ (4.11) 

If tuo >> -D, then we will have damped oscillatory motion. In the other extreme limit 
oi D » ujq, the motion becomes overdamped. While one solution of the second-order 
differential equation suffers strong damping, the other independent solution is neghgibly 
damped . The physical reason for this is that, under strong friction, any oscillator displaced 
from equilibrium and released from rest has only to acquire a small "terminal" velocity, so 
that friction balances the driving force. An oscillator starting from this "phase" of oscillation, 
will then almost freeze, and the associated energy in the oscillator will decrease negligibly. 
On the other hand, an oscillator mode with a large initial velocity will be significantly 
damped by the strong friction. Therefore, it is important to consider the ratio D/ujq to 
determine which limit applies for the nonlinear Alfven modes. 

We focus primarily on damping by photon viscosity. This is the most important source 
of viscosity, after e"'"e~ annihilation. Also, it is the dissipative process with the potential 
to damp the largest scales. (Smaller-scale damping by neutrinos is briefiy discussed in 
Section VII and Appendix C). The radiative viscosity coefficient is given by Eq. (p.9|), with 
g = g^ = 2, and the photon mean-free-path is 

, . X 1 9^ f T \-^ f n^h^ y^ , , , 

r = — ^ 9.5 X lO^^cm —^ x;\ 4.12 

^^ ' arneir) \0.25eV J VO-0125; " ^ ' 

Here, ay is the Thomson cross-section for electron-photon scattering, rif. is the electron 
number density, x^ the ionisation fraction, and Qb is the baryon density of the universe 
Pb, in units of the closure density. (For later convenience we define the quantity /;, = 
(r2b/0.0125/i~^), which measures the baryon density, in terms of the preferred value given 
by Walker et al. [|l^. This used to be the canonical value determined from nucleosynthesis 
constraints, although at present there is some debate on this issue |1^). In the early universe. 



the energy density of the baryon-photon fiuid is dominated by the photon energy, density 
p^; so, using p* = p*/3, p = p* /a^, tj* = o?rj, we have for the damping coefficient 

_ r,*e _ 3 4 P,{r)l,{T)e _ 1 
^ - (p* +p*) - « 15 (4p,(r)aV3) " 5^ ^'^^'^'' ^^>' ^^'^'^^ 

where we have used the fact that l^y^r) oc ri^^(T) oc a^i^r) and, as before, tq is the conformal 
time today. Hence, the damping-to-driving ratio is 

9 



D _ r/*fc2 ^ I ^^(^^);^(^) 



UJQ kVA{p*+p*) 5 Va 

^ 526.3 ^"^^^^"^^^ (4.14) 

-D-9 

We have defined the proper wavenumber kp{T) = {k/aij)), of a Fourier component and 
substituted for the Alfven velocity in terms of the field strength using Va ~ 3.8 x 10~^i?_9. 
For the diffusion approximation to be valid, we require kpl^ < 1; that is, we must consider 
only wavelengths larger than the mean-free-path. Nevertheless, one expects a large range of 
wavelengths for which modes will fall in the overdamped regime. 

In order to consider the evolution of modes of different wavelengths, one has first to look 
at the quantitative solution of Eq. ( [4.10| ) . For this, substituting 



fir) = exp (- 1 ^dr) W^(r) (4.15) 

into Eq. ( |4.10| ), the evolution of W is given by 

W + p{t)W = with p{T)=ul-j-{jy. (4.16) 



When ujq » D, we have p ~ cUq in Eq. ( |4.16| ) and then the solution is W = exp{±iujQT) . 
Therefore, in the oscillatory limit we have 

6o(r, z) = exp I- f ^-^dr] e±*-o-+^'=^, uJo» D. (4.17) 



In the opposite limit, D » ujq, we have to solve the oscillator with a time- dependent 
friction coefficient. One can obtain an approximate WKBJ solution. 



^(^) = -(Z^^^P 



± {-p)y\T)dT 



(4.18) 



This solution is valid as long as p(r) does not vary too rapidly. In the overdamped regime, 
this condition can be shown to be equivalent to neglecting D compared to D^. In the limit 
D « D'^, the two solutions are given by 
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fir) = Ao-^exp [-fD{r')dr') ; /(r) = i?o^exp [- f ^fr'^j (4.19) 

As advertised, in the overdamped limit, one solution (the Aq - mode) is strongly damped 
while the other solution (with Bq ^ 0) is weakly damped. 

For damping by photon viscosity, we have D oc a^ oc [t/tq)'^ in the radiation-dominated 
epoch. It is more useful to consider an alternate treatment to that of the WKBJ solution, for 
the rapidly varying, strongly overdamped regime. One notes that, as the damping increases 
with time to -D >> uq, f will tend to adjust itself so that the acceleration vanishes, so / ~ 0. 
For example, consider initially the case / > 0, / < and Df > uj^f. Then / > and so 
the magnitude of / decreases (while remaining negative) until we have Df = —uJqJ, when 
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the acceleration vanishes. On the other hand, if Df < uJqJ, then / < and the magnitude 
of / increases until we have / = 0. Subsequently as D keeps increasing, / can continue to 
adjust itself to maintain zero /. One can argue similarly for all other cases. Therefore, it 
seems plausible to consider an alternate approximation for the overdamped case, whereby, 
after the time when / first vanishes, / satisfies the equation 



Df + u^,f = 0; D»uJo. 



(4.20) 



We will refer to this approximation as the terminal-velocity approximation. The solution in 
the overdamped regime, under the terminal-velocity approximation, is simply given by 



/(^) = /(^T)exp 



iOn 



D(t' 



-dr' 



(4.21) 



Here, tt is the conformal time when the mode reaches the terminal- velocity regime, or when 
the acceleration, /, first vanishes. As we explained earlier, for an oscillator with an initial 
phase such that / is already large, this implies strong damping by the time the terminal- 
velocity regime is reached. On the other hand, for an oscillator which starts from rest, / 
will have to increase negligibly for the D » uq regime, so that / vanishes and Eq. (f4.21|) 
applies. 

We now move from the study of the non-linear Alfven modes to consider the damping 
of all the different MHD modes in their linearised limit. This has already been studied by 
JKO P). However, we shall do this using the formalism developed here to bring out the 
links with the non-linear situation. We will return to the further evolution of the non-linear 
Alfven mode in Section VI. 



V. DAMPING OF LINEARISED COMPRESSIBLE MHD WAVES 



Let us begin with the linearised MHD equations describing small perturbations to density 
pI = p* — Pq, pressure pi = p* — Pq, and magnetic field b* = B* — Bg. We have 



dr 



+ v.[(pS + p>]=o. 



(5.1) 



d 



„ , [V X b*] X B* 



VV + -V(V.v) 



(5.2) 



dh* 

'd7 



V X [v X B^ 



(5.3) 



V.b* = 0, pS = f, Pl = j- 



(5.4) 



Suppose the perturbation is described in terms of the perturbed comoving position (5x = 
^(x, r). The perturbed velocity v = {d^/dr) = ^. An integration of the perturbed continuity 
equation (|5J] ) and induction equation (^]3|) then gives 
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Substituting Eq. 
1 



Pi = -gPoV.e b* = V X [e X B*] . 
into the perturbed Euler equation, we get 



3?7* 



e = ^V(V.O + [V X (V X [^ X v^])] X v^ + 
3 4po 



v'^ + -v(v.o 



(5.5) 



(5.6) 



where we have defined va = Vaz. This hnear equation describes the evolution and damping 
of hnearised MHD modes in the expanding universe. One can look for plane wave solutions 
of the form ^ = ^/^(r) exp(ik.x). This leads to a replacement of V by ik and leads to an 
evolution equation for the amplitude ip'- 



^ 



(k.^)[, 



Vj] 



VAkz^Jz 



- Viki^/j - kzik.^p)Viz - 



3r]* 
4pS 



k'^ip + -kik.ip) 



(5.7) 



Here a subscript "z" denotes the 2;-component of the relevant quantity. We can now look at 
various types of solution to the above equation. 

First, consider the incompressible mode, with V..^ = 0, or k.^ = 0. In this case taking 
the dot product of Eq. (|5.7| ) with k, we also have ipz = (provided kz = kcosO ^ 0). For 
this mode, Eq. (|5.7|) reduces to 



3r)* 

V' + -^k^^ + vl cos^ ee^ = 0. 



(5.8) 



This is, as expected, is almost exactly the evolution equation for the Alfven mode encoun- 
tered and discussed in detail in the last section. The only difference is that Va is replaced 
by V^cos^, generalising the Alfven mode propagation to be in a general direction inclined 
at an angle 6 to the zero-order magnetic field. (Although this generalisation is at the cost 
of introducing the linear approximation). 

The evolution of the compressible modes, can be derived by taking the dot product of 
Eq. ( |5.7| ) with k and z. Defining A = k.ip/k, we have 



A + 



Vljk^A - Vie cos" e^z + \k^A 

Po 



0, 



(5.9) 



1 

— / 
3 



ijz + -k^ cos M + —k^z + —k^ cos 6 A = 0. 



4p5 



4p$ 



(5.10) 



Consider first the undamped limit with t]* = 0. In this case, looking for modes with A ex 
ipz oc e*"^"^, we can easily derive the dispersion relations 






cl + Vi 



±1 
2 



[cl + V^f-^cMcos'e 



1/2 



(5.11) 



Here, we have defined the sound speed in the relativistic limit Cs = l/v3. The plus sign in 
the above equation corresponds to the fast MHD mode while the negative sign corresponds 
to the slow mode. In the limit Va << Cs, which is generally applicable to our early universe 
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context, the dispersion relation for the fast mode becomes uo/k 
mode becomes oj /k ~ Va cos 6. 



C, while that of the slow 



The general solution of Eq. ( |5.9| ) and ( p.lO| ) when the damping terms are reinstated is 
quite complicated to analyze analytically, since it involves a fourth-order differential equation 
with time-dependent coefficients and the dispersion relation is a fourth-order poljTiomial. 
However, we can look at some simple special cases which illustrate the general behaviour. 

First, consider the case where k is parallel to Bq; then cos^ = kz/k = 1, A = ipz-, and 
the equation for A reduces to 



A + 'Lk'^A + elk'' A = 0. 
Po 



(5.12) 



This describes a damped sound wave, well studied in the literature in connection with the 
Silk damping of acoustic baryon-photon fluctuations. Therefore we only look at it briefly, 
to estimate the Silk damping scale. Because the sound-wave oscillation frequency is such 
that Us = kcs >> kVA, these modes do not become overdamped in general, and the damped 
oscillatory solutions of Section IV can be used to describe their evolution. Specifically, we 
have 



A(t) = exp 



Ds 



T 



-dr 



,,HlJsT 



UJ. 



»Ds 



(5.13) 



where Dg = k'^ij]* / p^). These modes get damped by a factor 



exp 



DJr) 



-dr 



exp 



^2 



where k^' 



2 
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l/fy \Aj Li 



(5.14) 



This agrees quite well with the Silk damping of sound waves in the radiation era, derived 
in more detailed treatments (cf. |T^, [0), in the appropriate limit. In fact, the more 
detailed derivation of Silk damping, using a Boltzmann treatment gives a damping factor 
exp[ 



-(^V^z),6o/)]> where 



'^D,bol 



l^dt 



15 J a^{t) 



[l + R)f-^ + 5Ry4 

(i + ny 



(5.15) 



Here, R = 3pb/4:p^ and /2 = 3/4 if the effects of polarisation are included and /2 = 1 
otherwise. Notice that, in the limit R << 1 applicable to epochs where the radiation density 
dominates baryon density, and /2 = 1, the damping scales exactly match, with k^ = ko^boi- 
The effects of non-zero R constitute at most about a 20 — 25% correction to the damping 
scale we derive. In the radiation-dominated epoch one has k]^^ = {A/AhY/''Ls{t)/a{t) ~ 
0.3L5(t)/a(t), where Lg{t) = (l^tY^'^ is the Silk scale. The largest scales which suffer 
appreciable damping are the modes with wavelengths (27rfc5^)) of order Ls [0. 
In the other extreme case, when k is perpendicular to Bq, cos 6' = k^/k - 
equation for A reduces to 



0, and the 



A + l-k'A+[- + Vi]k'A = Q. 
Po -J 



(5.16) 
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This describes a damped fast-magnetosonic wave. The real part of the oscillation frequency 
is uji = k{c^ + V^Y^'^ >> kVA, since in general Va << Cg. One can see that the damping 
of these modes is very similar to that of the sound waves (and in fact is exactly the same 
when we neglect Va compared to c^). Again, we expect modes with wavelengths less than 
the Silk scale to be significantly damped. 

Now we turn to the damping of modes with arbitrary direction k. In the undamped case, 
with rj* = 0, and when Va << Cg, we have seen that the fast mode has the same oscillation 
frequency as the sound wave and the slow mode as the Alfveh wave. This suggests an 
approximation to capture the damped counterparts of these modes in the limit of weak 
magnetic field. Let us write the time variation of A and ipz as A{t) cc ip^ oc exp {i J udr). 
Suppose u is dominated by its real part, and this is of order the undamped frequency of 
oscillation. We shall later check the consistency of this assumption. Consider first the 
damped counterpart of the slow mode, in the hmit Va « cj. = 1/3. For this mode, the 
ratio of the first two terms in Eq. ( |5.9| ), is A/{c^k'^A) ~ V^/cl « 1. Also, the ratio of the 
last term in (|0| ) to the second is ~ DA/{c'^k'^A) ~ VAkpl^/5 « 1. So, these two terms 
can be neglected when compared to the second term in Eq. (|5.9| ), so c^k'^A ^ V^A;^ cos Oipz- 
For the same reason, one can neglect the last term in Eq. ( |5.10|) compared to the second 
term in this equation. Substituting c'^k'^A ^ V^k"^ cos 9ipz in Eq. ( |5.10| ), we then have for 
the damped counterpart of the slow mode, 

3/7* 

i^z + -r^k'^^z + Vie cos^ e^z = 0. (5.17) 

We see that this is exactly the same equation as that obtained for the damped Alfven mode, 
analysed in the last section. Therefore, the slow modes will also be overdamped, and have 
one solution with negligible damping rate. Our original assumption that uj is dominated 
by its real part is valid for this solution, showing the self-consistency of our assumptions. 
Also, for this mode we have A ^ (Vl/cj.) cosOipz « i^z, for a general 9; that is, the mode 
is almost incompressible. (The strongly damped mode has to be analysed differently). 

Now consider the fast mode in a similar fashion. We have already derived the exact 
evolution for the special case when cos 6* = 0. Suppose cos 6* 7^ 0; assume, as before, that 
UJ is dominated by its real part, and this is of order the undamped frequency of oscillation. 
Then, substituting for the time dependence of A and ipz, and taking the real part of Eq. 
( p.l(j| ), we have ipz ~ A cos 6. So the ratio of the ipz term in Eq. ( |5.9D compared to the 
second term is ~ (V|/c^) cos^ 9 « 1. Neglecting the ipz term, and neglecting V^ compared 
to c^ in Eq. (|5.y|), results once again in the same equation, Eq. ( p.l2| ), for the fast wave as 
was found for the damped sound waves. 

In summary, the above analysis for linear MHD waves, in the limit of weak fields with 
Va << Cs, shows that the fast magnetosonic waves generally damp like sound waves, while 
there is one mode of the slow magnetosonic wave which behaves exactly as the Alfveh mode 
and gets overdamped. This also agrees with the conclusions reached by JKO. 

VI. THE FREE-STREAMING REGIME 

As the universe expands, the mean- free-path of the photon increases as a^, while the 
proper length of any perturbed region increases as a. So the photon mean-free-path can 
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eventually become larger than the proper wavelength of a given mode. When this happens 
for any given mode, we will say that the mode has entered the free-streaming regime. Modes 
with progressively larger wavelengths enter the free-streaming regime up to a proper wave- 
length ~ l'y{Td) ~ 10^^ cm (see Eq. (|4.12| ) ), or a comoving wavelength of ~ 3Mpc, at the 
epoch of decoupling. After (re) combination of electrons and nuclei into atoms, l^ increases 
to a value larger than the present Hubble radius, and all modes enter the free-streaming 
regime. (We will consider the pre- and post-recombination epochs separately below.) 

When photons start to free stream on a given scale of perturbation, the tight-coupling 
diffusion approximation no longer provides a valid description of the evolution of the per- 
turbed photon-baryon fluid on that scale. One has to integrate the Boltzmann equation 
for the photons together with the MHD equations for the baryon-magnetic field system. A 
simpler approximate method of examining the evolution of such modes in the linear regime 
is to treat the radiation as isotropic and homogeneous, and only consider its frictional damp- 
ing force on the fluid. (The radiative flux could have also contributed to the force on the 
baryons; however, for modes with wavelengths smaller than l^, this flux is negligible since the 
associated compressible motions have suffered strong Silk damping at earlier epochs; when 
the wavelength was larger than l^). The drag force on the baryon fluid per unit volume due 
to the radiation energy density p^, is given by 

4 

Fd = --r2eO-TP7V. (6.1) 

Since, typically, less than one electron-photon scattering occurs within a wavelength, the 
pressure and inertia contributed by the radiation can be neglected when considering the 
evolution of such modes. The Euler equation for the baryonic component then becomes 

-^ + H{t)Y + V. Vv = Vpb + — J X B - -V0 - -^Uearv. (6.2) 

ot apb Pb ct 3pb 

Here, pb is the baryon density, pb the fluid pressure, and H{t) = {da/dt)/a is the Hubble 
parameter. We have also included the gravitational force, (l/a)V0, due to any perturbation 
in the density. Note that we have written this equation in the unstarred conformal frame, 
(with the magnetic field defined in the "Lab" frame; see Appendix A), since conformal 
transformation to flat spacetime is no longer a useful tool in the matter-dominated era. We 
have also transformed the time co-ordinate, from conformal time, back to "proper time" 
dt = adr. 

It should be pointed out that the dramatic drop in the pressure, by a factor of order the 
very small baryon to photon ratio ~ 10~^, when a mode enters the free-streaming regime, 
has important consequences. First, in the absence of radiation pressure, the effect of mag- 
netic pressure (if it greatly exceeds the fluid pressure) is to convert what was initially an 
incompressible Alfveh mode into a compressible mode (see below). Second, the effective 
baryonic Jeans mass decreases dramatically and compressible modes can become gravita- 
tionally unstable. Thus, we have to retain the gravitational force term in the above equation. 
The magnetic pressure will also play a dominant role, providing pressure support against 
gravity on sufficiently small scales. 

The evolution of modes which enter the free-streaming regime depends on the strength 
of the magnetic fields, in particular whether the magnetic pressure pb is greater or smaller 
than the fiuid pressure, pb- For the magnetic pressure, we have 
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Pb = ^{l + zy^Ax 10-^B% (^) dyn/cm\ (6.3) 

While the fluid pressure is given by 

p, = 2nefcT ^ 1.1 X 10-1° (^^) hdyn/cm\ (6.4) 

where we have assumed that the fluid temperature is locked to the radiation temperature, 
and that the gas is an electron-proton gas. By taking the ratio of the two pressures, one can 
see that magnetic pressure dominates the fluid pressure, (i.e. ps » Pb ioi B » Bcrit ~ 
5 X 10^11 Gauss). For magnetic fields smaller than Bcru, the fluid pressure dominates. 

Consider first the case where the field B is much smaller than Bcru- In this case the 
motions can be assumed incompressible. The Alfven modes which enter the free-streaming 
regime, remain Alfvenic. Following the ideas of Section IV, we look again at non-linear 
Alfven modes with B = (Bq + b)/a^, where Bq = Bqz, with Bq = constant, b = nbo{z,t) 
and V = nvo{z,t), with n perpendicular to z. Recall that |b| is not necessarily small 
compared to |Bo|. We assume pb to be uniform (but not independent of t), use the Euler 
equation ( |6.2D and the induction equation (PT7|), change to conformal time r, and look for 
solutions in the form bQ{z,t) = f{t)e^^^, following the same procedure as in Section IV. (For 
the rotational Alfven- type mode, the gradient terms in (|6.2| ) do not contribute). We obtain 
an equation for the evolution of /(r) 



d'f . \ rr . r<\df 



dr' 
where 



+ 



aH + D :^+uj^=0, (6.5) 

dr 



u;o = kVA{^Y/'; D = n,aTa{^). (6.6) 

Here, V'^ = 3-BQ/(167rp^) = constant, is the same as defined in terms of the starred variables 
in Section IV. 

The evolution of this non-linear Alfven mode depends once again on the relative strengths 
of the damping and driving terms. The ratio of the viscous damping to expansion damping 
is given by 

D ^ i4:py/3pb)neaTa ^ ^Dh_ ^^ ^ ,^ ^. 

aH aH 3pb l^ 

since the Hubble radius Dh = H~^ » l^, so one can neglect the damping due to Hubble 
expansion. Also, the ratio of the viscous damping to the driving terms in the oscillator 
equation ( |6.5| ) is 

5 ^ iAp,/3pb)n,aTa _ 3 .^ iq^(P^)U2 1 .g o^ 

uo kVA{Ap,/3pbY/^ ■ V/ k,it)l,it)B^,- ^ • ^ 

When a given mode enters the free-streaming limit we will have kp{t)l^{t) ~ 1. So, for the 
field strengths -B_9 < (i?cr«t/10^^G) << 1 that we are considering, all the Alfven modes are 
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strongly overdamped. As the universe expands, the product kp(t)l^{t) oc a^ increases, and 
at late times any given mode enters the damped oscillatory regime. One can again apply 
the terminal-velocity approximation of Section IV in the overdamped regime once df/dr 
has adjusted itself to the "zero acceleration" solution. Then, / is given by 




/(r) = Cexp -/ -^dr']; D » u^. (6.9) 



At sufficiently late times, such that D /ujq << 1, the damped oscillatory solutions are ap- 
propriate, so 

/ = exp ( - / ^^dr] e^^^«", uo^ » D. (6.10) 



Now consider the case B » Bcru, when magnetic pressure dominates over the fluid 
pressure. For p^ >> ph, incompressibility is no longer a good assumption, and we necessarily 
will also have gravitationally unstable, magneto-acoustic modes. In general, the pressure- 
gradient term (dominated by the magnetic fleld) and the magnetic tension are of similar 
order. The non-uniform magnetic flelds (associated with what were initially Alfven-type 
modes), can seed compressible motions and their associated density fluctuations. These 
density fluctuations, can grow by gravitational instability into non-linear structures. Before 
recombination, there is still signiflcant drag due to the free-streaming photons. But once 
atoms form, the density of free electrons decreases and the photon mean-free-path becomes 
larger than the Hubble radius. In the free-streaming regime, the evolution equations of 
these gravitationally unstable MHD modes are the same before and after recombination. 
We consider a unifled analysis of this evolution in the next section. 

VII. THE POST-RECOMBINATION REGIME 

We assume that the perturbations in density and velocity are small enough so that non- 
linear terms in the perturbed density and velocity can be neglected. In the Euler equation 
( |6.2D , one can neglect the non-linear term, v.Vv and take the density pb to be that of 
the unperturbed FRW background density. This equation has to be supplemented by the 
continuity equation for the perturbed fluid density, the Poisson equation for the potential, 
and the induction equation (|3.7|) . These equations are given explicitly in Appendix B. In 
the Poisson equation, we take account of the possibility that there may be other forms of 
coUisionless dark matter by writing 

VV = 47rGa^5pr = 47rGa^ [pbSb + pJc] ■ (7.1) 

Here, Spx is the total perturbed density (due to both fluid plus dark matter), Sb is the 
fractional perturbation in the fluid, while pc and Sc describe the dark matter density and 
its fractional perturbation, respectively. We shall adopt the equation of state pb = 2nekT = 
Pb{2kT/mp) = pbcf, where rup is proton mass. 

In general, the magnetic fleld will be non-uniform when the mode enters the free- 
streaming regime. We note that, when the background fluid pressure was large,the non- 
uniform magnetic fleld may have originally been part of an Alfven-type incompressible mode. 
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However, as we mentioned earlier, once this mode enters the free-streaming regime, there is 
a dramatic fall in the fluid pressure, by a factor of order the very small baryon to photon 
ratio ~ 10~^. As a result, the pressure of the non-uniform magnetic field, associated with 
what might well have been an Alfven-type mode, can no longer be ignored, especially if the 
field exceeds Bcru- (Only a perfectly circularly-polarised Alfven wave has uniform magnetic 
pressure). This non-uniform field associated with what started off as an Alfven-type mode, 
will now also induce gravitationally unstable, compressible motions. 

In treating the resulting evolution, it is usual to assume (cf. Wassermann [|1^] and 
Peebles jT^) that perturbations to the Lorentz force, due to that the perturbed velocity, are 
subdominant with respect to the zeroth-order contribution of the Lorentz force itself. So one 
takes B = Bo(x)a^, which solves the induction equation (p.7|) , if v is neglected. Of course, 
this approximation will break down once significant peculiar velocities have been developed, 
as will always happen on sufficiently small scales, or at sufficiently late times, for any given 
magnetic field. For galactic scales, it turns out that the distortions to the magnetic field 
will become important only at late times, even for B_g ~ 1. So the above assumption of 
retaining only the zeroth-order contribution to the Lorentz force is expected to be reasonable. 
(The equations governing the more general case are derived in Appendix B, and we hope 
to return to an analysis of this full system elsewhere). Making these assumptions, standard 
linear perturbation analysis (cf. ||16i, |T^ , [IH]), leads to the evolution equation for 5^, 



dH, 



' + 



4p 

2H H -UeCTTa 



-^ - clV% = AnGa^ [ptSb + pM + -j5o(x) (7.2) 



where the source term 5*0 is given by 

_ V. [Bq X (V X Bq)] 

Here, pb{to) is the fiuid density at the present time, to. If we assume the dark matter to be 
cold, one can also derive a similar equation for its fractional perturbed density 6c- One finds 

^ + 2H^ = AnGa' [p,6, + p,5e] . (7.4) 

Prior to recombination, as pointed out in the previous section (Eq. (|6.7|) ), the viscous 
damping dominates damping by the Hubble expansion. For Bq > Bent, the fiuid pressure 
term can be neglected. Let us also assume that the magnetic field was the only source of 
initial density perturbations. Then, the compressible modes start their free-streaming evolu- 
tion with negligible initial (5^ since modes smaller than the Silk scale have been significantly 
damped and modes on larger scales have a negligible source of pressure perturbations in the 
radiation era because ~ 3V^/c^ ~ 3 x lO^^i^.g. Thus, Eq. ( [7.2| ) can be solved, under the 
terminal velocity approximation, simply equating viscous damping and the Lorentz force. 
We obtain 

Here, Ib is a typical co-moving coherence scale over which the field varies. Note that the 
RHS of this equation is constant in time and so the density contrast 6b increases linearly 



in this epoch. Hence, at the time of recombination, tr, the induced baryonic perturbation 
is 6b = ■Jditr — tf), where tf is the time when the scale k~^{t) ~ iBait), becomes smaller 
than the photon mean-free-path. For a flat universe, dominated by dark matter, the total 
fractional density contrast in the matter is, 6 = {ptSb + Pc6c)/{Pb + Pc) ~ ^b^b- At the time 
of recombination, we then have 

S{U) ^ 3.8 X 10-'B%h-' y^^ (1 - tj/U). (7.6) 

On galactic scales with Ib ~ IMpc, we have tf/tj. ~ 0.4, and then 5{tr) ~ 2.1 x 10~^h~^B'^g. 
This turns out to be small compared to the 6 induced in the post-recombination regime (see 
below) . 

Next, consider the post-recombination evolution. The mean-free-path of the photon now 
increases rapidly to a value exceeding the Hubble radius and viscous damping becomes 
subdominant compared to expansion damping. Thus, we can neglect the viscous damping 
term. Also, for B > Bcrit, we can neglect the fluid pressure term. Now suppose the baryons 
contribute a fraction fs to the matter density while the cold dark matter contributes a 
fraction 1 — fs- Then, multiplying Eq. ( [7.2|) by /b, ( [r.4|) by 1 — /b, and adding the resulting 
equations, we get for the total density contrast 6 = {pbSb + PcSc)/{Pb + Pc), 

Here, pm is the total matter density. Let us assume that at recombination 5 ~ and 
{85/ dt) ~ 0; that is, initially there are negligible fluctuations in density and velocity diver- 
gence. Note that this is valid for scales much larger than galactic scales. For scales with 
Ib ~ Mpc, it turns out that, the post recombination evolution induces a 5 much larger 
than that given in Eq. ( [7.6|) , within an expansion time. So including this initial 5{tr) gives 
negligible corrections. The particular solution of ( |7.7| ) for a fiat matter-dominated universe 
is given by 



m = ^fstiso 






(7.8) 



This implies that the magnetic field induces a present-day fractional density contrast 6o = 
(5(to), with 

So = l^fBtlSoil + ..) - ^/b (^^M^y (1 + ,^) (7.9) 

where (1 + Zr) = (t/t^)^/^ and V^6(to) is the Alfven velocity with respect to the baryons 
given by 



VAbito) Bf) 



^47rpb(to)c2 



1.9 X lO-^B-g/^r'^' (7.10) 



Adopting h = 1/2, 1 + Zj. = 1100, and a fiat universe, with to = 2/(3iJo)) Ib = 0.05, we 
have 
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*-2.96By^)- (7.11) 

We see therefore that a magnetic field with Bq ~ lO^^G is needed to impact significantly on 
galaxy formation. Such fields will also induce rotational perturbations and give significant 



angular momentum to protogalaxies [0. 

iFTom Eq. ( [7l8|) and ( [7.11|) , it would seem that for sufficiently small, Ib = h, say, one 
can have a 6b ~ 1, even close to recombination. However, the calculation leading to Eq. 
( [7.8|) would break down on such small scales because the field distortions induced by the 
motions, which we have neglected, will become dominant. The resulting magnetic pressures 
will oppose gravity when the Alfven crossing time becomes of order the dynamical time; that 
is, for proper lengths a{t)lB < hit) ~ VAb(t)t, where t is the relevant dynamical time (the age 
of the universe) and lj{t) is the magnetic Jeans length. Noting that VAb{t) = VAbito)a~^^'^{t) 
and t = toa^^"^ for a fiat universe, the comoving magnetic Jeans length is given by 

Aj = ^ ~ VAb{to)to ~ 3.8 X 10~^B^,h~'f^'/^Mpc. (7.12) 



Our treatment of how inhomogeneous magnetic fields induce structure formation is valid 
only for Ib » \j- On smaller scales one has to solve the full set of equations outlined 
in Appendix B. We expect strong magneto-sonic waves to be induced by such small scale 
inhomogeneities in a sufficiently strong magnetic field. These may suffer strong dissipation 
and so input energy into the IGM. We hope to return to this issue elsewhere. 

VIII. DISCUSSION 

We have studied the evolution and damping of inhomogeneous magnetic fields in various 
regimes. It is of interest to synthesise our results and discuss how a given spectrum of 
magnetic inhomogeneities evolves. We make a few general points and then describe the fate 
of magnetic inhomogeneities after they enter the Hubble radius. 

• As we noted in Section IV, for magnetic field strengths B^g < 1, the pressure pertur- 
bations are negligible. The fast compressible waves induced by these perturbations, 
have a phase velocity of order the relativistic sound speed, but negligible amplitudes 
with (5fe ~ Va/^'s ~ 3y|/c^. The motions induced by the field can then be treated as 
incompressible, at least until they enter the free-streaming regime. These slow residual 
motions occur no faster than the Alfven timescale, and pressure can constantly read- 
just on the fast sound-crossing time scale to preserve the incompressibility condition. 
In this case, Alfven modes, both non-linear and linear, and the incompressible limit of 
the slow mode, that we studied in previous sections, are indeed the most relevant. 

• It is important to note that linearisation about a constant background field is not a 
good approximation when following the evolution of magnetic fields which are inho- 
mogeneous with roughly similar power on a multitude of scales. This was one of our 
motivations for concentrating on the non-linear Alfven mode, and studying its evolu- 
tion and damping through various epochs. Although this analysis employed special 
exact solutions, the amplitude of the tangled component of the magnetic field could be 
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taken to be arbitrarily large compared to the amplitude of the large-scale field. Also, 
its spatial configuration can be arbitrarily specified, by the free function bo{TQ, z). So, 
one expects the behaviour of this mode to refiect, at least qualitatively, the behaviour 
of general incompressible motions driven by magnetic-field inhomogeneities. 

A comoving scale Ib-, which enters the Hubble radius in the radiation-dominated epoch, 
does so at an epoch tg, specified by /Ba(te) = Stg- Suppose we define this epoch by the 
radiation temperature T{t) = To/a(t), where Tq is the present-day microwave background 
radiation temperature, then we have fir 



nte)=n = 6^eV{^)-\ (8.1) 

In models which produce the field during an infiationary epoch, the initial condition for the 
non-linear Alfven wave at the time of horizon entry could be taken to be that the fiuid is 
at rest but the field is tangled. The Lorentz force due to this tangled field will then start 
pushing on the fiuid, when the scale of the tangle becomes smaller than the Hubble radius. 
On the other hand, if the fields are produced in an early-universe phase transition, they 
could be associated with large initial velocities. 

• In the cosmological context, it should be kept in mind that there is only a finite time 
for the Alfven wave to develop and induce motions in the fiuid, if they were initially 
absent. For example, on a comoving scale of k~^, by the time of recombination, at 
r = Tr, (or at a temperature T = T^), the Alfven wave would have oscillated at most 
by a phase angle of 

X = kV,(r.. - r.) < kV.n ~ 1.8 x 10-B_,(^)-' (^)"' . (8.2) 

Here, Te is the conformal time when a mode enters the Hubble radius, and we have 
expressed x iii terms of the inverse of the wavenumber, k~^. Thus, only small-scale 
magnetic inhomogeneities, with scales k~^ < /^ ~ O.lB^gMpc, will have had time 
to oscillate by more than it/2 in phase. Modes on scales k~^ » Is, which started 
with zero initial peculiar velocity, and oscillated by a phase x << 1 when damping 
is ignored, cannot damp their tangles, even if damping is included. This situation 
applies to magnetic fields tangled on galactic scales, with k~^ ~ Mpc. Thus galactic- 
scale magnetic inhomogeneities, with B^g < 1, do not get damped by photon viscosity, 
simply because the incompressible wavelike motions they induce, oscillate negligibly 
before recombination. 

• As the photon mean-free-path grows, the fast compressible motions, induced by the 
magnetic field or due to existing "adiabatic" density fiuctuations, are damped on scales 
smaller than the Silk scale, that is on scales k~^ less than about 0.3Ls ~ 0.3(/^(t)t)^'^. 
The comoving Silk damping scale, L^ = Ls/a, at any time t (or temperature T) in 
the radiation era, is given by 

L%{T) ^ 8.5 X 10^^ [-^^^) K"'"cm. (8.3) 
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After matter domination, when the scale factor a(t) oc t^/^, the co-moving Silk scale is 
given by 

L%{T) ^ 6 X 10^^ \^l^ h-'/'f-'/'cm. (8.4) 

By contrast, the evolution of non-linear Alfven wave modes depends on the ratio D/ujq. 
Non-linear Alfven wave modes which enter the Hubble radius when the temperature 
of the universe T > 2.5eV, do so when D/ujq < 1, and the mode is then initially in the 
damped oscillatory regime. But the photon mean-free-path grows as a^ while wave- 
lengths grow as a, and soon tangles on some scale are in the overdamped regime, with 
D/ujq > 1. For a given scale, k~^, this happens at a time tqd, when the temperature 
of the universe drops to T < Tqd, where 

Too = lO.SeV X B-Ti^J-'^'fV'^"- (8-5) 

So, tangles on smaller scales not only enter the Hubble radius at an earlier epoch, but 
are overdamped at an earlier epoch. 

While the Alfven mode is in the damped oscillatory regime, its evolution can be de- 
scribed by Eq. ( [4.17|) , and so all the modes damp by a factor of order exp(— J'^{D/2)) ~ 
exp[— (3/4)(4/Cp(r)L|(r)/45)]. This is almost the same as the Silk damping factor, for 
the usual baryon-photon sound waves. Hence, all modes which are smaller than the 
Silk scale will get significantly damped by photon viscosity before the non-linear Alfven 
mode enters the overdamped regime. The largest comoving scale, say k~^ = k~g]., which 
gets damped while the mode is in the damped oscillatory regime, can be estimated by 
equating [koscLs(ToD,kosc)]/'^5 = 1. This gives 

k-^ ^ 4.4 X IQ-^ B'tjbMpc. (8.6) 

In the damped oscillatory regime, a mode which starts initially from rest (/ = 0), 
can be described by the solution /(r) ~ /o cos[x(t — Tg)] exp[— J{D/2)]. Even modes 
which damp negligibly in this above regime, will oscillate by a phase of order 

xiroD) - kVAToD ^ 4.14 X lO~^B%\^r'/'fl/' (8.7) 

when the mode enters the overdamped regime. Modes with xiToo) > 1 will acquire a 
large / by the time they enter the overdamped regime and will be strongly damped. 
This will happen for any mode with k^^ = kQ^^ < 1.7 x lO^^B'^gf^ Mpc. 

It turns out that, for the Alfven-type modes with k~^ > S/c^d, which enter the Hubble 
radius with zero initial velocity, there is negligible further damping while the mode 
is in the overdamped regime. When these modes enter the overdamped regime, with 
ujq/D = 1, we find that they do not have sufficient velocity (/) for friction to be 
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important. With increasing time, the velocity grows, D increases, and ujq/D de- 
creases. Eventually, the mode enters the terminal-velocity regime, at say a time tt, 
where friction balances the Lorentz force (see Section IV). This happens roughly when 
tan(x(rr)) ~ {ojq/ D{tt)) « 1 and /(t^) ~ /o- In the terminal velocity regime, we 
have from Eq. (|4.21|) , that 



/(^) = /(^T)exp 



uo, 



TT 



ZJ* 



/(rr) exp 



-0.9 X 10" 



,Tt-T 
'0.25e1/ 



)Bl 



(8.8) 



For all modes with Ib > 51od, one finds that there is no significant damping. 

When the photon mean-free-path increases above the wavelength of a given mode, one 
enters the free-streaming regime. As we discussed in Section VI, the further evolution 
of the Alfven-type mode depends on the whether magnetic pressure at this time domi- 
nates the fluid pressure, or vice versa. This is determined by the ratio B/Bcru- For the 
case of B << Bcru, one can again treat the evolution as incompressible. The resulting 
Alfven modes are already in the overdamped regime when they start to free-stream 
and their evolution is governed by Eq. (|6.9|), 



/(^) = /(^/)exp 



oo' 



^/ 



-s,t 



/(r;)exp-Vi 



tf 



kl{t)l^{t)dt 



(8.9) 



Hence, /(r) = /(r/) exp(— A;^/A;j^), where the free-streaming damping scale k^} is 
given by 



klj = V I '^ 



l^{t)dt 



a\t) 



(8.10) 



Modes with a scale for the magnetic field k~^ < kj^, get damped significantly during 
the free streaming evolution. We see that the damping in this regime is similar to Silk 
damping, except that the usual Silk damping integral within the exponential (cf. Eq. 
( p. 141 ) ) is multiplied by an extra factor of (15/2)V| << 1. For the linearised modes 
where the wave vector makes an angle 6 to the zero-order field, one has to replace 
Va by V^cos^. After recombination, the viscous damping is subdominant, compared 
to expansion damping (since l^ exceeds the Hubble radius), and so can be neglected. 
So the largest scale to be damped is found by evaluating kj^ at the recombination 
redshift. Assuming that the universe is matter dominated at recombination, we get 
kj^ ^ {3/5Y^'^VALg{tr). Hence, the damping scale is of order the Alfven velocity 
times the Silk scale. The largest wavelength mode to be damped, say L^ = 2TrkJ^{tr) 
is given by. 



- - ^'(!'"^^^ 



1.1 X 10'''B^,f,-'/'h-'/^cm. 



(8.11) 



For B > Bcrit, we noted in Section VI that the evolution becomes compressible, and 
gravitationally unstable for scales larger than the magnetic Jeans length, Xj. On such 
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scales, we showed that a field with _B_9 ~ 1, is needed to produce a density pertur- 
bation large enough to significantly affect galaxy formation. We did not treat the 
evolution on scales smaller than the magnetic Jeans length in any detail, although 
the governing equations are given in Appendix B (cf. Eq. (B8)). The solution of 
these equations, as we noted in Section VII, is complicated by the presence of an inho- 
mogeneous zero-order magnetic field. Nevertheless, we expect that fast compressible 
motions on scales smaller than Xj will drive oscillations close to the baryonic Alfven 
frequency, and will be initially overdamped by the action of photon viscosity, in the 
pre-recombination era. The damping scale for such motions will then be similar to 
kj^, as deduced above. For modes similar to the slow magneto-sonic wave, we expect 
that for B » Bcrit, the phase velocity will be roughly equal to Cb, the baryon sound 
speed, and the corresponding damping scale will be kj^ with Va replaced by Cfe. These 
expectations are borne out by the linearised calculations of JKO, although they ignore 
the inhomogeneous nature of the zero-order magnetic field. (For a general tangled zero 
order field, the counterpart of the slow wave may not be easy to excite). Clearly, more 
detailed computations are needed to get the exact damping scales, in this case. 
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• One can consider the damping due to neutrino viscosity in the early universe ||23 
an exactly analogous manner to the Silk damping effects treated in detail above. We 
briefly describe below some of the consequences of damping due to neutrinos. These 
damping effects are largest around the time of neutrino decoupling, at a temperature 
of T = T,^ ~ IMeV. During this epoch, the number density of weakly interacting 
particles is riw ~ T^ (cf. [|l9l, |l^). The cross-section for interaction with neutrinos, 
is typically the weak interaction cross-section aw ~ G^^T^, where Gp is the Fermi 
constant. The neutrino mean- free-path is then given by 

l^ = — i — ^lAxl0^^cm(-f—) . (8.12) 

awnw \MeVj ^ ' 

Note that l^ becomes comparable to the Hubble radius, H~^ = 2t ^ 4.4 x 
10^^ cm{T/MeV)~'^, at T ~ IMeV. The comoving damping scale for fast modes, 
say Lg'^, due to neutrino viscosity, can be derived in exactly analogous manner to 
the Silk damping scale, derived in Section V. We get Lg'^ = [t{T^)l^{T^)Y/'^ /a{T^) ~ 
ly{Ty) / a{Ty) ^ 0.9 X 10^°cm. As with the damping by photon viscosity, the maxi- 
mum damping of Alfven (or the slow) modes occur when they enter the free-streaming 
regime. We analyse the damping of the non-linear Alfven wave in the neutrino free- 
streaming regime, in Appendix C We show there that damping in this regime is similar 
to Silk damping except that the usual Silk damping integral within the exponential 
(cf. Eq. ( |5.14|) ) is multiplied by an extra factor of (15/2)V|y << 1, and l^ is replaced 
by l„. Here, Vau is the Alfven velocity defined in terms of the conserved, neutrino 
energy density py{Ty)a^{Ty). Since l^ oc T^^ oc t^^^ in the radiation-dominated epoch, 
we can estimate the largest wavelength Alfven mode that is appreciably damped, by 
neutrino viscosity, in the free streaming regime. This comoving wavelength is given 
by Li, ~ 27r(2/5)VV^,Lg'^ ~ WB_,cm. 

We can generalise these results given any initial spectrum of magnetic inhomogeneities, by 
replacing V^ by a suitably averaged scale-dependent V^(A;~^). 
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IX. CONCLUSIONS 

We have considered the evolution and viscous damping of cosmic magnetic fields in the 
early universe in detail. Using the fact that the fiuid, electromagnetic, and shear viscous 
energy-momentum tensors are all conformally invariant, we showed in Section II that the 
MHD equations in the FRW universe, including viscous effects, can be transformed into their 
special-relativistic counterpart when the metric perturbations from inhomogeneous motions 
are small. This enabled us to transform known non-linear Alfven-wave solutions, from flat 
spacetime, into the expanding FRW universe. 

We considered in detail the evolution and damping of these modes in various regimes. 
First, on galactic scales or larger, the Alfven mode oscillates negligibly before recombination, 
for magnetic fleld strengths, B^g < 1 (or a present day magnetic fleld B < 10~^ Gauss). So 
there is then no question of strong "Silk" damping of these modes, due to photon viscosity, 
as occurs for compressional baryon-photon oscillations. Furthermore, Alfven waves with 
small enough wavelength, which can oscillate appreciably before recombination, become 
overdamped. In this case, the longest wavelength which suffers appreciable damping by 
photon viscosity, has a scale k~^ ~ VaL^, where Lg is the usual comoving Silk scale. Since 
the Alfven speed is Va ~ 3.8 x 10~'^B-g << 1, only comoving wavelengths smaller than 
10^^i?_9 cm suffer appreciable damping. We also briefly considered analogous results for 
damping of very small scales by neutrino viscosity following neutrino decoupling at t ~ Is. 

After recombination, the fluid pressure drops enormously, roughly by the baryon to pho- 
ton ratio. The Lorentz force due to a tangled magnetic fleld associated with an initially 
Alfven-like mode before recombination can seed gravitationally unstable compressional per- 
turbations after recombination, provided the fleld is strong enough and tangled enough on 
scales larger than the magnetic Jeans length, Xj. We examined the post-recombination evo- 
lution of scales larger than Xj in section VII, including the effect of a passive dark matter 
component, and showed that magnetic flelds with 5_9 ~ 1 are needed to impact signifl- 
cantly on galaxy formation. The evolution equations for perturbations with scales smaller 
than Xj are derived in Appendix B, but their solution is much more complicated and will 
be examined elsewhere. 

Magnetic flelds with B_g ~ 1, may be constrained by observations of quasar rotation 
measures [^. The fluid velocities induced by the tangled fleld, oscillating as an Alfven 
mode in the pre-recombination era, can also produce small angular scale anisotropies in the 



microwave background, through the Doppler effect |^2[- This may provide another constraint 
on such flelds. The dissipation of magnetic flelds due to neutrino or photon viscosity will 
also leave an imprint on the spectrum of neutrinos and photons respectively. The neutrino 
spectrum could be probed using nucleosynthesis and there already exist strong limits on the 
spectral distortions of the microwave background. Scenarios of galaxy formation that appeal 
to strong enough magnetic flelds of order 10~^ Gauss are potentially testable. We hope to 
return to a further consideration of some of these observational issues in future work. 
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APPENDIX A: MAXWELL'S EQUATION IN THE "LAB" FRAME 

In the main text we defined the electric field E* = [E*^ , E*^ , E*^) and magnetic field 
B* = (-B*\ B*^, B*^) in the starred metric by 

This can be written in more compact way in term of a four-vector electric field E* and 
magnetic field B* as 

B; = le;,,xV*''F*''\ E; = F;,V*^, (A2) 

where V*'-^ = [1, 0, 0, 0] is the four-velocity of fundamental observers at rest in the starred 
metric frame. We have also used the Levi-Civita tensor e*^^;^ = ^/—g*A^up\, with A^upx, the 
totally antisymmetric symbol such that ^0123 = 1 and ±1 for any even or odd permutations 
of (0,1,2,3). Note that the four-vectors B* and E* have purely spatial components and 
E* = £■**, B* = B*^ for the spatial components of the field 4- vector. First, let us transform 
the electric and magnetic four-vectors to the unstarred (FRW) conformal frame, and denote 
the resulting electric and magnetic four-vectors by E'^ and B'^, respectively. Making use of 
the conformal transformation properties. 



we get 



.*,.^^-4^,u. ,*^^^ = ^4^^^^^. y- = fi-V^ (A3) 



B'^ = nB; = ^B;; e'^ = qe; = ^e;. (A4) 



We now make a coordinate transformation to the FRW proper "Lab" co-ordinates (t, r) 
defined by dt = adr, dr = ad:>c. Denote the co-ordinate-transformed, "Lab" electric and 
magnetic four- vectors by E^^ and i?^, respectively. We have 

b, = -b'^ = \b;; e, = -e'^ = \e;. (A5) 

Note that these 4- vectors are also purely spatial. Similarly, we can define the current density 
in the Lab frame by J^ = aJ^ = aVt'^J*^ = a~^J*^. Now, denote the spatial components 
-B^, E^ and J^ by the spatial 3- vectors B, E and J respectively. Then, in terms of these 
spatial vectors, we have 

B* E* J* ,. , 

B = — ; E=— ; J = - A6 

a^ a^ a-^ 

The four-current density will also have a time component, the charge density Pq = Jl = 
{J*^/a^). Using the Maxwell equations, ( p.4|) , in the starred metric, one can write the 
Maxwell equations in terms of these " Lab" fields. We have 
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V X (a^B) = 47ra^J + 



dr 



V.B = 



(A7) 



V X (a'E) 



dr 



V.{a^E) = Ana^pg. 



(A8) 



The non-relativistic limit of Ohm's law in terms of the "Lab" fields is simply 

E + vxB = -. (A9) 

a 

In the ideal limit, the "Lab" magnetic field then satisfies the induction equation 

V X (a'B) 



d{a'B) _ld{a'B) ^1^^ 



dt 



a dr 



(AlO) 



So, when v = 0, we have B oc a~^, a result which is intuitively expected for the "Lab" 
magnetic field due to flux freezing in the expanding universe. It is also of interest to express 
the Euler equation ( p.3|) for the fluid in the radiation-dominated era, in terms of the "Lab" 
fields. We obtain 



]_d_ 



(p + p)a'^v 



fv.V) V 

+ ^^ [(P + P)v] + -V. [(p + p)v] 

1 n 

= — Vp + JxB + 4 



V^v + -V(V.v) 



(All) 



APPENDIX B: GENERALISED LINEAR PERTURBATION THEORY IN THE 
FREE-STREAMING AND POST-RECOMBINATION REGIMES 



Begin with the linearised evolution equations for the baryonic fluid, including the effects 
of the magnetic field and gravity. We have 



-^ + H{t)^r = — -Vpi + -J X B - -V0 - ^n^arv 
ot aph pb a 3pb 



(Bl) 



dt 



-V.v = 0. 



(B2) 



V^ 



AtxGo?5pt = 47rGa^ [pbh + Pc^c 



(B3) 



Here pi = cfpbSb is the perturbed pressure. We take B = [Bo(x) + b(x)]/a^, with |b| << 
|Bo|. The linearised induction equation, (|A10D, then becomes 



dh 

'dt 



-V X [v X B 







V.b = 0. 



(B4) 
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Suppose the perturbation is described in terms of the perturbed position of the baryonic 
component by 5x{, = ^(x, r). The perturbed velocity is v = a{d^/dt). An integration of the 
perturbed continuity equation (|E^ ) and induction equation ( P^ ) then gives 

(B5) 



-V.e b = V X [e X Bo] . 



Similarly, one can define the perturbed comoving position of the cold dark matter component 
by 5xc = ^d^iT)- An integration of the perturbed continuity equation for this component 
gives 5c = —'V.C,c- The Poisson equation ( jB^ ) for the potential can then be integrated once 
to give 



- -V0 = AnGa [pb^ + p^^ 
a 

The perturbed pressure-gradient term can be written as 

- — vp, = ^v(v.o. 

apb a 



(B6) 



(B7) 



Using Eq. 



( |B6| ) and ( |B7| ), the perturbed Euler equation becomes 



+ 



2H + ^n^ar 
^Pb 



iV X (V X [e X Bo])] X Bo 



(47rpba^)a^ 
(V X Bo) X [(V X [e X Bo)] 



{A7Tpha^)a^ 



(B8) 



This linear equation describes the gravitationally unstable evolution, and damping of lin- 
earised MHD modes in the free-streaming and post-recombination regime. However, one 
cannot perform a simple Fourier analysis of this equation, since Bo is also a function of x. 
In the main text we have simplified the equation by neglecting the last two terms. While 
this approximation is likely to be valid for large-scale modes at early times, it will break 
down once the distortion of the magnetic field due to the motions become significant. We 
hope to return to the study of this equation elsewhere. 



APPENDIX C: 



VISCOUS DAMPING DUE TO NEUTRINOS IN THE 
FREE-STREAMING REGIME 



Modes whose wavelengths become longer than the neutrino mean-free-path, enter the 
neutrino free-streaming regime. The evolution and damping of non-linear Alfven wave so- 
lutions can be examined in this regime, analogous to the case of photon-free streaming, 
treated in section VI. First the viscous force due to coupling of the relativistic plasma and 
the neutrinos is given by |2^ 



'd = -^nwawpuv = ^Pi^Y' 



(CI) 
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Here, p^ is the energy density of neutrinos, and we have replaced {nwcrw)~^ by the neutrino 
mean-free-path ly. We can also use the form of the Euler and induction equations, (as 
derived in Appendix A); except that the viscous force is as given by Eq. ( Pl| ) and the 
fluid inertia/pressure does not include the neutrino contributions. Following the ideas of 
Section IV and VI, we look again at non-linear Alfven modes with B = (Bq + b)/a^, where 
Bo = -Bqz, with Bq = constant, b = ng{T)e^'^^ and v = nVg{T)e'^'^^, with n perpendicular to 
z. Recall that |b| is not necessarily smaW compared to |Bo|. Also, since the universe is still 
radiation dominated at the epochs near neutrino decoupling, with p » B'^/{87t), we can 
treat the motions as nearly incompressible. The Euler equation gives 

Apna^dv Bodh 4 v 

where pn is the fluid density, excluding the neutrino contribution. The induction equation 
gives 

c9b dv 

These can be combined to give an equation for the evolution of (^(r) 

-f^ + Dy— + u^g = 0, (C4) 

(Jj I Hi I 

where 

u, = kVAR, D. = f(^)ocr-^ (C5) 

l-u ipR 

Here, V\ji = 3-Bo/(167rpija^) = constant, is the Alfven velocity defined with the inertia 
contributed by p^. 

The evolution of this non-linear Alfven mode depends once again on the relative strengths 
of the damping and driving terms. When modes enter the neutrino free-streaming regime, 
the Alfven waves are again strongly overdamped, with Dy/uju >> 1, for B_q < 1. One can 
again apply the terminal-velocity approximation of Section IV in the overdamped regime, 
once dg/dr has adjusted itself to the "zero acceleration" solution. Then, g is given by 

g{T) = g{n) exp {- f k^VlJ-^dr] , (C6) 

\ Jr,, a 



where 



B, 



2 



^^'^ - 4vr(4plv3) ^^^^ 

is the Alfven velocity defined in terms of the conserved neutrino energy density pyO^. We 
see that the damping of these Alfven modes by neutrino viscosity is similar to Silk damp- 
ing except that the usual Silk damping integral within the exponential (cf. Eq. ( |5.14 )) is 



multiplied by an extra factor of (15/2)V2y << 1, and l^ is replaced by the neutrino mean 
free-path l^. 



29 



*Presently at: Max Plank Institute for Astrophysics, Karl Scliwarzschild Strasse 1, 85748 
Garching, Germany. On leave from National Centre for Radio Astrophysics, TIFR Pune, 
Poona University Campus, Ganeshkhind, Pune 411 007, India. 



30 



REFERENCES 

[1] Moffat, H. K., Magnetic Field Generation in Electrically Conducting Fluids, Cambridge 
University Press, Cambridge, (1978); Krause, F. & Radler, K.-H., Mean-Field Magne- 
tohydrodynamics and Dynamo Theory, Pergamon Press, Oxford, (1980); Ruzmaikin, A. 
A., Shukurov, A. M. & Sokoloff, D. D., Magnetic Fields of Galaxies, Kluwer, Dordrecht 
(1988); Mestel, L., Stellar Magnetism, Oxford UP, Oxford, (1998); Beck, R., Branden- 
burg, A., Moss, D., Sliukurov, A. & Sokoloff, D., Ann. Rev. Astron. Astropliys., 34, 155 
(1996). 

[2] Vishniac, E.T., Ap. J. 446, 724 (1995); Rees, M. J., Quart. Jl. Roy. Astron. Soc. 28, 
197 (1987); Rees, M.J., In Gosmical Magnetism, ed. Lynden-Bell, D., Kluwer, London, 
(1994) pl55; Subramanian, K., Bull. Astr. Soc. Ind., 23, 481 (1995); Subramanian, 
K., Narasimha, D. & Cliitre, S.M., Mon. Not. R. astron. Soc. 271, L15 (1994), and 
references therein. 

[3] Hogan, C. J., Phys. Rev. Lett., 51, 1488 (1983); Turner, M. S. & Widrow, L. M., Phys. 
Rev. D, 30, 2743 (1988); Vachaspati, T., Phys. Lett. B, 265, 258 (1991). Ratra, B.,. 
Ap.J Lett., 391, LI (1992); Gasperini, M., Giovannini, M. & Veneziano, G., Phys. 
Rev. Lett. 75, 3796 (1995); Lemoine, D., & Lemoine, M., Phys. Rev. D, 52, 1955 
(1995); O. Tornkvist, |hep-ph/97075T3| ; K. Enqvist, |astro-ph/9707300| ; Giovaninni, M., 



and Shaposhnikov, M.E., |hep-ph/9710234| ; Joyce, M. and Shaposhnikov, M.E., Phys. 
Rev. Lett., 79, 1193 (1997). 

[4] Cattaneo, F. & Vainshtein, S. I., Ap.J., 376, L21 (1991); Vainshtein, S. & Rosner, 
R., Ap.J., 376, 199 (1991); Kulsrud, R.M. & Anderson, S.W., Ap.J., 396, 606 (1992); 
Brandenburg, A., In Lectures on Solar and Planetary Dynamos, ed. Proctor, M. R. E. 
& Gilbert, A. D., Cambridge University Press, (1994), pll7; Field, G. B., The Physics 
of the Interstellar Medium and Intergalactic Medium, ed. Ferrara et al., ASP Conf. 
Ser., 80, (1996), pi; Blackman, E., Phys. Rev. Lett., 77, 2694 (1996); Chandran, B., 
Princeton University Observatory preprint (1996) POPe-659; Cattaneo, F. & Hughes, 
D. W., Phys. Rev E, 54, 4532 (1996); Subramanian, K., Mon. Not. R. astron. Soc. in 
press; |astro-ph/9707280| . 

[5] Kulsrud, R., lAU Symp. 140: Galactic and Extragalactic Magnetic Fields, Reidel, Dor- 
drecht, (1990), p527, and references therein. 

[6] J.D. Barrow, P.G. Ferreira, and J. Silk, Phys. Rev. Lett. 78, 3610 (1997); J.D. Barrow, 
Phys. Rev. D 55, 7451 (1997). 

[7] Silk, J., Ap. J., 151, 431 (1968). 

[8] Jedamzik, K., Katalinic, V., and Olinto, A., (1996), |astro-ph/9606080| v2 

[9] Weinberg, S., Gravitation and Cosmology, Wiley, New York, (1972); G.F.R. Ellis, in 
General Relativity and Cosmology, ed. R.K. Sachs, Varenna Lectures in Physics, Aca- 
demic, New York, (1971). 
[10] Brandenburg, A., Enqvist, K., Olesen, P., Phys. Rev. D, 54, 1291 (1996). 
[11] Liang, E. P. T., Ap. J., 211, 361 (1977). 

[12] Parker, E. N., Cosmic Magnetic Fields, Clarendon Press, Oxford (1979). 
[13] Landau, L. & Lifshitz, E. M., Electrodynamics of Continuous Media, Pergamon Press, 

Oxford (1987). 
[14] Walker, T. P., Steigman, G., Schramm, D. N., Olive, K. A. & Kang, H.-S., Ap.J., 376, 
51 (1991). 

31 



[15] Steigman, G., Hata, N. & Felten, J. E., |astro-ph/9708016] and references therein. 

[16] Wasserman, L, Ap.J., 224, 337 (1978) 

[17] Peebles, P. J. E., Large Scale Structure of the Universe, Princeton University Press, 

Princeton N. J., (1980). 
[18] Padmanabhan, T., Structure Formation in the Universe, Cambridge University Press, 

Cambridge (1993). 
[19] Kolb, E. W. & Turner, M. S., The Early Universe , Addison- Wesley, New York, (1990) 

[20] Kaiser, N.„ Mon. Not. R. astron. Soc, 202, 1169 (1983). 

[21] Kronberg, P. P., Rep. Prog. Pliys., 57, 325 (1994) 

[22] Subramanian, K. & Barrow, J. D., 1998 (in preparation). 

[23] Neutrino viscosity was first considered as a means of damping anisotropics and inliomo- 
geneities by C.W. Misner, Nature 214, 40 (1967) and Pliys. Rev. Lett. 19, 533 (1967). 
See also, J.D. Barrow, Nature 267, 117 (1977). The damping of very small scale irreg- 
ularities by GUT scale transport processes is discussed by J.D. Barrow, Mon. Not. R. 
astron. Soc. 199, 45P (1982). 



32 



